First Results from the Cryogenic Dark Matter Search in the Soudan Underground Lab 
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We report the first results from a search for weakly interacting massive particles (WIMPs) in the 
Cryogenic Dark Matter Search (CDMS) experiment at the Soudan Underground Laboratory. Four 
Ge and two Si detectors were operated for 52.6 live days, providing 19.4 kg-d of Ge net exposure after 
cuts for recoil energies between 10-100 keV. A blind analysis was performed using only calibration 
data to define the energy threshold and selection criteria for nuclear-recoil candidates. Using the 
standard dark-matter halo and nuclear-physics WIMP model, these data set the world's lowest 
exclusion limits on the coherent WIMP-nucleon scalar cross-section for all WIMP masses above 
15 GeV/c^, ruling out a significant range of neutralino supersymmetric models. The minimum of 
this limit curve at the 90% C.L. is 4 x 10"*^ cm^ at a WIMP mass of 60 GeV/c^. 



PACS numbers; 95.35.-|-d, 95.30.Cq, 14.80.Ly 

There is a compelling scientific case that nonlumi- 
nous, nonbaryonic, weakly interacting massive particles 
(WIMPs) 0,1111 may constitute most of the matter 
in the universe |^. Supersymmetry provides a nat- 
ural WIMP candidate in the form of the lightest su- 
perpartner, which must be stable if R-parity is con- 
served IE 11 S II H- The WIMPs are expected to be 
in a roughly isothermal halo whose gravitational poten- 
tial well contains the visible portion of our galaxy. These 
WIMPs would interact elastically with nuclei, generating 
a recoil energy of a few tens of keV, at a rate smaller than 
-1 event kg-^ d'^ H .S. li IS H E3 ■ 

The Cryogenic Dark Matter Search ( CDMS) collabo- 
ration is operating a new apparatus [llj to search for 
WIMPs in the Soudan Underground Laboratory. The 
CDMS Soudan experiment, also called CDMS II, uses 
a set of Ge (each 250 g) and Si (each 100 g) ZIP (Z- 
dependent Ionization and Phonon) detectors , cooled 
to temperatures < 50 mK and surrounded by substan- 



tial shielding deep underground to reduce backgrounds 
from radioactivity and cosmic-ray interactions. Simulta- 
neous measurement of ionization and athermal phonon 
signals in the Ge and Si detectors allows excellent re- 
jection of the remaining gamma and beta backgrounds. 
These background particles scatter off electrons in the 
detectors, while WIMPs (and neutrons) scatter off nu- 
clei. The ZIP detectors allow discrimination between 
electron and nuclear recoils through two effects. First, 
for a given energy, recoiling electrons are more ionizing 
than recoiling nuclei, resulting in a higher ratio of ioniza- 
tion to phonon signal, called "ionization yield." Second, 
the athermal phonon signals due to nuclear recoils have 
longer rise times and occur later than those due to elec- 
tron recoils. For recoils within a few /im of a detector's 
surface (primarily from low-energy electrons) , the charge 
collection is incomplete , making discrimination based 
on ionization yield less effective. But these events can be 
effectively rejected by phonon timing cuts because they 



2 



have, on the average, even faster phonon signals than 
those from bulk electron recoils [l^ Il5| . These effects are 
in qualitative agreement with our understanding of the 
complex phonon and semiconductor physics involved |l6l| . 

The detectors are surrounded by an average of 0.5 cm 
of copper, 22.5 cm of lead, and 50 cm of polyethylene, 
which reduce backgrounds from external photons and 
neutrons. A 5-cm-thick scintillator muon veto enclos- 
ing the shielding identifies charged particles (and some 
neutral particles) that pass through it. An overburden 
of 780 m of rock, or 2090 meters water equivalent (mwe) , 
reduces the surface muon flux by a factor of 5 x 10"*. 

All materials surrounding the detectors have been 
screened to minimize radioactive decays which could pro- 
duce neutrons. Neutrons resulting from radioactive de- 
cays outside the shield are moderated sufficiently to pro- 
duce recoil energies below our detector threshold. Neu- 
trons produced in the shield by high-energy cosmic-ray 
muons are tagged by the veto scintillator with an ef- 
ficiency >99%. The dominant unvetoed neutron back- 
ground is expected to arise from neutrons produced by 
cosmic-ray muon interactions in the walls of the cav- 
ern. Events due to neutrons can be distinguished in part 
from ones due to WIMPs because neutrons often scat- 
ter in more than one detector and interact at about the 
same rate in Si and Ge. By contrast, WIMPs would not 
multiple-scatter, and coherent scalar WIMP interactions 
would occur ^6x more often in Ge than in Si detectors. 

We report here the analysis of the first CDMS Ge 
WIMP-search data taken at Soudan during the period 
October 11, 2003 through January 11, 2004 [uj. After 
excluding time for calibrations, cryogen transfers, main- 
tenance, and periods of increased noise, we obtained 
52.6 live days with the four Ge and two Si detectors 
of "Tower 1" (six close-stacked ZIP detectors labeled as 
Zl(Ge), Z2(Ge), Z3(Ge), Z4(Si), Z5(Ge) and Z6(Si) from 
top to bottom). Tower 1 was operated previously in an 
identical configuration at the Stanford Underground Fa- 
cility (SUF), at a depth of 17 mwe 

Energy calibrations were performed repeatedly during 
the run using a ^'^^Ba gamma source with distinctive, 
penetrating lines at 356 keV and 384 keV. The excellent 
agreement between data and Monte Carlo simulations 
and the observation of the 10.4 keV Ga line from neutron 
activation of Ge indicated that the energy calibration was 
accurate and stable to within a few percent. Observation 
of the predicted energy spectrum from a ^^^Cf source 
confirmed the energy scale for nuclear recoils [ll| . 

The trigger rate on phonon signals for the WIMP- 
search data runs was ^0.1 Hz, with a recoil-energy 
threshold of ~2 keV (~4 keV for Zl). The muon veto 
signals were recorded in a time history buffer for each 
detector trigger. The summed veto rate above threshold 
was ~600 Hz, mainly due to ambient gammas. At this 
rate, we reject about 3% of events with accidental veto 
coincidences in the 50 ^s before a detector trigger. Data- 
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FIG. 1: Ionization yield versus recoil energy for calibration 
data with a ^""^Cf gamma and neutron source for detectors Z2, 
Z3 and Z5 in Tower 1 showing the ±2(7 gamma band (solid 
curves) and the ±2(j nuclear-recoil band (dashed curves) for 
Z5, the detector with the worst noise of these three. Events 
with ionization yield < 0.75 (grey) are shown only if they 
pass the phonon-timing cuts. The vertical line is the 10 keV 
analysis threshold for these three detectors. 



quality cuts reject the ^5% of events that show any sign 
of higher pre-trigger noise or possible pile-up. 

We performed a blind analysis, in which the nuclear- 
recoil region for the WIMP-search data was not inspected 
until all cuts and analysis thresholds were defined using 
in situ gamma and neutron calibrations (see Fig. A 
combination of ionization-yield and phonon-timing cuts 
rejects virtually all calibration electron recoils while ac- 
cepting most of the nuclear recoils. The phonon timing 
cuts are based on both the phonon rise time and the 
phonon start time relative to the ionization signal (see 
Fig. 121. We required recoil energy between 10-100 keV 
for all Ge detectors except Zl, whose larger noise required 
an analysis threshold of 20 keV in order to ensure compa- 
rable rejection. We rejected events with some ionization 
in a detector's annular "guard" electrode, which covers 
15% of the detector's volume. Figure |21 shows a conser- 
vative estimate of the combined efficiency of all cuts on 
a WIMP signal. The cuts yield a spectrum-averaged ef- 
fective exposure of 19.4 kg-days between 10-100 keV for 
a 60 GeV/c^ WIMP. 

Table ^ lists the observed rates of unvetoed events in 
the WIMP-search data, with ionization yield either in 
the ±2(7 gamma band ("gammas") or below this band 
(mostly surface electron recoils). Analysis shows that 
about half the surface electron recoils with interactions 
in only a single detector ( "singles" ) were due to beta de- 
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FIG. 2: Phonon start time versus ionization yield for Ba 
gamma-calibration events (diamonds) and ^^^Cf neutron- 
calibration events (dots) in the energy range 20-40 keV in 
detector Z5 in Tower 1. Lines indicate typical timing and 
ionization-yield cuts, resulting in high nuclear-recoil efficiency 
and a low rate of misidentified surface events. 
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FIG. 3: Efficiency of the combined cuts as a function of recoil 
energy, both for the blind analysis (solid) and for the second, 
non-blind analysis (dashed). The step at 20 keV is due to 
Zl's 20-keV analysis threshold. 



cays of contaminants on surfaces, while the other half 
were from gamma rays. Gamma rates are ~50% higher 
at Soudan than they were at SUF, consistent with the 
higher Rn levels at Soudan and the absence of a 1-cm- 
thick ancient-Pb liner which surrounded the detectors at 
SUF. Total surface-event rates at Soudan are also some- 
what higher than at SUF, consistent with the increased 
component due to gammas. 

We computed the number of electron-recoil events 
expected to be misidentified as nuclear recoils in the 
WIMP-search data based on the ^"^^Ba calibration sets 
used to determine the timing cuts. Factoring in system- 
atic errors, we estimated 0.4 ±0.3 misidentified events in 
Zl and a total of 0.3 ±0.2 in the other Ge detectors. As a 
check, we applied the same cuts to a different set of ^'^^Ba 
calibrations, containing 1.5 times as many surface events 
as in the WIMP-search data. One event (at 50 keV in Zl) 



TABLE L Unvetoed gamma and surface-electron-recoil rates 
between 15-45 keV in Tower 1 at Soudan. 





Gammas [#/kg/day] 


Surface [^^/day] 


ZIP 


(total) 


(singles) 


(total) 


(singles) 


Zl(Ge) 
Z2(Ge) 
Z3(Ge) 
Z5(Ge) 


85.6 ± 3.4 
79.4 ± 3.1 
89.3 ± 3.3 
105.7 ± 3.6 


37.6 ± 2.3 

19.7 ± 1.6 
19.9 ± 1.5 
35.9 ± 2.1 


1.56 ± 0.23 
1.05 ± 0.18 
1.11 ± 0.18 
1.82 ± 0.24 


0.90 ± 0.18 
0.18 ± 0.08 
0.15 ± 0.07 
0.65 ± 0.14 
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FIG. 4: Ionization yield versus recoil energy for WIMP-search 
data from Z2 (triangle), Z3, and Z5 (+) in Tower 1, using the 
same yield-dependent cuts and showing the same curves as in 
Fig-0 Above an ionization yield of 0.75, the events from all 
three detectors are drawn as identical points in order to show 
the 10.4 keV Ga line from neutron activation of Ge. 



passed all cuts, in agreement with the previous estimate. 

Monte Carlo simulations predict 0.05 ± 0.02 neutrons 
(mostly unvetoed) produced from muon interactions out- 
side the shielding, including uncertainties on the neutron 
production rate. The simulations predict ~1.9 (veto- 
coincident) neutrons produced inside the shielding for 
the WIMP-search data. No veto-coincident nuclear-recoil 
candidates were observed in the WIMP-search data. 

This blind analysis of the first Soudan CDMS II 
WIMP-search data set revealed no nuclear-recoil events 
in 52.6 kg-d raw exposure in our Ge detectors. Figure^ 
displays the ionization yield of WIMP-search events in 
Z2, Z3, and Z5 which passed the same cuts applied to cal- 
ibration data in Fig.^ As shown in Fig.[Sl these data to- 
gether with corresponding data for Zl set an upper limit 
on the WIMP-nucleon cross-section of 4 x 10~^^ cm^ at 
the 90% C.L. at a WIMP mass of 60 GeV/c^ for coherent 
scalar interactions and a standard WIMP halo. 

After unblinding the nuclear-recoil region, we found 
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regions under other assumptions see |22|)- Our new lim- 
its significantly constrain supersymmetric models under 
some theoretical frameworks that place weak constraints 
on symmetry-breaking parameters {e.g. 0.00I)- 
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FIG. 5: New limit on the WIMP-nucleon scalar cross section 
from CDMS II at Soudan with no candidate events in 19.4 
kg-d effective Ge exposure (solid curve). Parameter space 
above the curve is excluded at the 90% C.L. These limits 
constrain supersymmetry models, for example Q (dark grey) 
and (hght grey). The DAMA (1-4) 3cr signal region 0| is 
shown as a closed contour. Also shown are limits from CDMS 
at SUP [13 (dots), EDELWEISS 19] (x's), and the second, 
non-blind analysis of CDMS II at Soudan with 1 nuclear-recoil 
candidate event (dashes). All curves ^3 are normalized fol- 
lowing 10], using the Helm spin-independent form-factor, 
scaling, WIMP characteristic velocity vo — 220 km s~^, mean 
Earth velocity ve = 232 km s~^, and p = 0.3 GeV c~^ cm~"^. 



that our pulse-fitting algorithm designed to handle satu- 
rated pulses had been inadvertently used to analyze most 
of the unsaturated pulses in the WIMP-search data. This 
algorithm gives slightly worse energy resolution than the 
intended algorithm. The limit in Fig.|31based on the blind 
analysis (solid line) correctly accounts for this effect. We 
have also performed a second, non-blind analysis, using 
the intended pulse-fitting algorithm and the same blind 
cuts, resulting in a 5% higher WIMP detection efficiency. 
This analysis resulted in one nuclear-recoil candidate (at 
64 keV in Z5), consistent with the expected surface-event 
misidentification quoted above. Figure |5l includes the op- 
timum interval \21] limit based on this second unbiased, 
but non-blind, analysis (dashed line). 

At 60 GeV/c^, these limits are a factor of four below 
the best previous limits set by EDELWEISS jl^, and 
a factor of eight better than our limit with the same 
Tower 1 at SUE 0] . These data confirm that events de- 
tected by CDMS at SUE and those detected by EDEL- 
WEISS were not a WIMP signal. Under the assumptions 
of a standard halo model, our new limits are clearly in- 
compatible with the DAMA (1-4) signal region 18] if it 
is due to coherent scalar WIMP interactions (for DAMA 
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